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Abstract 



In a recent article by Dickey, J. O., Marcus, S.L. and O. de Viron, 2011, the authors 
show evidences for correlations in the multi-decadal variability of Earth's rotation 
rate [i.e., length of day (LOD)], the angular momentum of the core (CAM), and nat- 
ural surface air temperature (SAT). Previous investigators have already found that 
', the LOD fluctuations are largely attributed to core-mantle interactions and that 
the SAT appears strongly anti-correlated with the decadal LOD. As the above au- 
thors note, the cause of this common variability needs to be further investigated and 
studied. In fact, "since temperature cannot affect the CAM or LOD to a sufficient 
^h *| extent, the results favor either a direct effect of Earth's core-generated magnetic field 
Q H i (e.g., through the modulation of charged-particle fluxes, which may impact cloud 

formation) or a more indirect effect of some other core process on the climate-or yet 
another process that affects both" . The main aim of the present research note is to 
further support the above results, using a proper continuous wavelet analysis, and 
. to evidence some early detected time-period relations existing between LOD and 
Q\ ! SAT data. 

o 



1 Introduction 



We recall here some brief introductive considerations of the argument under study, as well 
reported in the reference article by Dickey, J. O., Marcus, S.L. and O. de Viron, 2011. The 
length of day (i.e., the time needed by Earth to make one full rotation) fluctuates around 
its mean 24-h value, over a broad range of periods. Earth rotation is observed with a very 
high accuracy from precise geodetic techniques such as very-long-baseline interferometry, 
which allows for the determination of the length of day (LOD) at the 0.02-ms level (less 
than 1 cm at the equator). Earth rotation data collected before 1975 are the results of 
traditional optical observations; thus uncertainties during this period are significantly larger 
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than for current determinations. The observed fluctuation of the LOD can be interpreted, at 
seasonal time scales, in terms of angular momentum exchange between the solid Earth and 
the superficial fluid layers (i.e. the atmosphere and the ocean), without the need to invoke 
influences from the core. For the atmosphere, variations in global angular momentum are 
dominated by the effect of the winds (the wind term), with the moment-of- inertia changes 
(the pressure term) being an order of magnitude smaller (Hide and Dickey 1991; Rosen 
1993). The ocean contributions, being yet smaller than the pressure term, are significant in 
the closure of the angular momentum balance (Marcus et al. 1998). On the other hand, at 
decadal and longer time scales (multi-decadal) , the main source of rotational variation is 
the interaction between the mantle and core, as substantiated by the significant correlation 
between the low-degree zonal components of the magnetic field and the LOD (Jault and 
Le Mouel 1991). The decadal LOD variations (e.g., ~ 4 ms around 1900) are too large in 
amplitude to be explained by the atmosphere (~ 1 ms in amplitude). The long-term global 
mean Earth surface air temperature (SAT) is significantly anti-correlated with decadal and 
longer LOD (e.g., Lambeck and Cazenave 1976). Variability at periods of 60-80 yr has 
been well established in LOD (e.g., Jault and Le Mouel 1991; Roberts et al. 2007) and in 
the core angular momentum (CAM), using both observational data (Zatman and Bloxham 
1997; Dickey and de Viron 2009) and theoretical studies (Mound and Buffett 2007). 

It is well known that there are no direct observations of the fluid core motions. However, 
geomagnetic fields, generated by the dynamo inside the fluid core, can be used to infer 
the motion field and its velocity by utilizing several basic assumptions and the induction 
equation. In fact, if the solid Earth-core system is assumed to be isolated, the changes of 
the solid-Earth angular momentum, and thus its rotation speed and the associated LOD, 
can be deduced from changes of the CAM. The main result of the above authors is that 
LOD, CAM and corrected (i.e. for anthropogenic effects) global mean air temperature show 
strong evidence of multi-decadal variability, with a period in the broad range: 60-80- yr. Os- 
cillations in global temperatures with periods in the 65-70- yr range were originally reported 
by Schlesinger and Ramankutty (1994). Subsequent observational studies and simulations 
with coupled atmosphere-ocean models have found similar multi-decadal climatic modes, 
typically originating in the North Atlantic Ocean; however, the excitation source or sources 
of these oscillations have not been unambiguously identified (Knight 2009). The authors 
work suggests that the same core processes that are known to affect Earth's rotation and 
magnetic field (Roberts et al. 2007) may also contribute to the excitation of such modes, 
possibly through geomagnetic modulation on regional as well as global scales (Usoskin et 
al. 2008). For more details see the full article by Dickey, J. O., Marcus, S.L. and O. de 
Viron, 2011. 



2 Results 



Following the final comment of the authors: "Further work remains to be done, especially 
in linking common modes of variability between Earth's subsystems and better describing 
the physical connections between them", here we used, as starting point, the curves of Fig. 
1 in the original article (see the reported Figure 1), considering the -LOD behavior and the 
corrected temperature series from GISS beginning in 1880 and ending 2002.9. 
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Fig. 1. The GISS observed and corrected temperature series (black and red lines, respectively) 
beginning in 1880 (cf. Hansen and Lebedeff, 1987). The LOD and CAM are given in units of 
equivalent milliseconds, with sign reversed so that positive anomalies correspond to increased 
rotational speed of the solid Earth, and the temperature series are given in units of 0.1 C. (From: 
Dickey, J. O., Marcus, S.L. and O. de Viron, 2011). 

The authors main results are: "...the detection of a broadband variability centered at 78 
yr (common variability ranges from 67 to 86 yr from SSA method). Oscillations in global 
temperatures with periods in the 65-70 yr are well known. Our work suggests that the 
same core processes that are known to affect Earth's rotation and magnetic field may also 
contribute to the excitation of such modes, possibly through geomagnetic modulation of 
near-Earth charged particle fluxes that may influence cloud nucleation processes, and hence 
the planetary albedo, on regional as well as global scales." 

The main goal of the following wavelet analysis is to confirm the above global correlation 
between -LOD and Corrected Temperature (i.e. avoiding the anthropogenic component) 
at multi-decadal scales and furthermore to better investigate time-frequency correlations 
between the series in order to strictly correlate the processes involved, observing their modes 
temporal behaviors. For mathematical details on this wavelet analysis, see: G. Ranucci and 
S. Sello, 2007. 

The following Figure 2 shows the wavelet analysis performed on the first series: -LOD data 
Note that the y-axis is shown in logarithmic scale in order to better see all the frequency 
range of interest involved (here selected from: 1 yr to: 121 yr). Black contour lines are 
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WAVELET ANALYSIS 




1.5E+0C 
1.2E+0C 
1.0E+0C 
8.4E-01 
6.9E-01 
5.7E-01 
4.7E-01 
3.9E-0I 
3.2E-01 
2.6E-0I 
2.2E-01 
1.8E-01 
1.5E-01 
1.2E-01 
1.0E-01 



J I I I I I l I I I I I L 



j i j_ 



880 1900 1920 1940 1960 1980 2000 

Time (year) 



Fig. 2. Wavelet analysis of -LOD data 
boundaries of confidence regions at 95%. 

Main results: we confirm a dominant and time-persistent frequency band in the periodicity 
range: 50-92 yrs, with a central peak at 65 yr, clearly visible up to f 900. Another persistent 
frequency band is near 34 yr. Further, other less persistent frequencies are: 33.8 - 13.8 - 
8.41 - 3.7 - 3.0 and 1.7 yr. Some of these frequencies are well time localized and may be the 
object of further future investigations. 



The following Figure 3 shows the wavelet analysis performed on the second series: Corrected 
T data. Note that the scales of the wavelet are the same as Figure 1 but the contour power 
levels here are rescaled to the different intensity energy level of the series. Black contour 
lines are boundaries of confidence regions at 95%. 

Main results: we confirm again a dominant and time-persistent frequency band in the range: 
60-120 yrs, with a central peak at 78.25 yr, clearly visible in the full time range considered. 
Another persistent frequency band is near 30.8 yr. Further, other less persistent frequencies 
are now: 19.8 - 9.0 - 4.7 yr. Some of these frequencies are again well time localized but with 
interesting differences, both in relative intensity and/or time-frequency position from those 
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WAVELET ANALYSIS 
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Fig. 3. Wavelet analysis of Corrected T data 

of the -LOD case. These may be the object of further future investigations. A qualitative 
analysis of the main common periodicities detected both in -LOD and Temp series, leads to 
the preliminary conclusion that there is a general shift towards longer periodicities (from 
-LOD to Temp data) together with a phase shift of few years, probably related to the 
time response of atmosphere to changes of LOD. The excitation of similar modes shifted 
in time seems to exclude some external processes that affect Earth's core and climate 
simultaneously. In fact, the authors already note "an approximately 8-yr lag between the 
(negative) LOD and the corrected temperature; this lag agrees with the 8-yr lag between 
changes in Earth's rotational speed and surface geomagnetic field perturbations found by 
Roberts et al. (2007)". 



The following Figure 4 shows the wavelet analysis performed on the third series: -CAM 
data Note that the scales of the wavelet are the same as Figure 1. Black contour lines are 
boundaries of confidence regions at 95%. 

Main results: we confirm again a dominant and time-persistent frequency band in the long 
multi-decadal range. Here, due to the smoother and shorter CAM data, we were able to 
resolve, as maximum period, only the value: 55.7 yr, clearly visible in the full time range 



5 




WAVELET ANALYSIS 
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Fig. 4. Wavelet analysis of CAM data 

considered. Other significant periodicities are: 33.1 - 22.2 - 13.2 - 8.16 and 2.6 yr. Some of 
these frequencies are again well time localized but with visible differences, mainly corre- 
sponding to different time positions from those of the -LOD case. However, there are clear 
common periodicities, as expected, with -LOD data. These localized periodicities may be 
the object of further future investigations. 



3 Conclusions 



The above wavelet analysis results support the idea, stated in the article by Jean O. Dickey, 
Steven L. Marcus and Oliver de Viron, that some sources of the (natural) global temper- 
atures oscillations could be identified with the core processes that affect Earth's rotation 
and magnetic field. In particular, the excitation of similar modes shifted in time seems to 
exclude the possibility of external (e.g. solar) processes that affects the core and climate 
simultaneously. It seems more likely that some core processes, that induce mode excitations 
of rotation and magnetic fields, drive global air temperatures oscillations too. Through the 
multi-resolution wavelet analysis we can see that not only multi-decadal oscillations but 
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also some high energetic decadal and sub-decadal well time-localized periodicities appear 
somewhat correlated and thus we need further investigations in order to give a correct in- 
terpretation. The time-frequency approach, available e.g. with a proper wavelet analysis, 
could help further detailed investigations and interpretations on this interesting suggested 
relation. 
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